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Abstract 

°. 

We use anticommuting variables to study probability distributions of random variables, 
■ that are solutions of Langevin's equation. We show that the probability density always 

enjoys "worldpoint supersymmetry" . The partition function, however, may not. We find 
that the domain of integration can acquire a boundary, that implies that the auxiliary field 
has a non-zero expectation value, signalling spontaneous supersymmetry breaking. This is 
due to the presence of "fermionic" zeromodes, whose contribution cannot be cancelled by a 
surface term. This we prove by an explicit calculation of the regularized partition function, 
as well as by computing the moments of the auxiliary field and checking whether they 
Q-r satisfy the identities implied by Wick's theorem. Nevertheless, supersymmetry manifests 

itself in the identities that are satisfied by the moments of the scalar, whose expressions we 
can calculate, for all values of the coupling constant. We also provide some quantitative 
estimates concerning the visibility of supersymmetry breaking effects in the identities for 
the moments and remark that the shape of the distribution of the auxiliary field can 
influence quite strongly how easy it would be to mask them, since the expectation value 
^ ' of the auxiliary field doesn't coincide with its typical value. 
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1 Introduction 



Let us consider a stochastic process, of commuting random variables, x(t), that satisfies the 
Langevin equation 

Here t](t) is a Gaussian stochastic process, i.e. its correlation functions satisfy the relations 



(V(r)) = 

(v(n)v(T2) ■ "flfon)) = (^(1)^(2)) (^(3)^(4)) ■ ■ ■ (?77r(2n-l)^7r(2n)) 



(v{n)v{T2)) = S(n -r 2 ) ^ 



where the sum is over all permutations, tt, of the index values 1, 2, 3, ... , 2n. 

We shall assume that U(x(t)) is an ultra-local functional of x(r), i.e. doesn't contain any 
derivatives with respect to r and discuss in our conclusions what happens when U(x(t)) is, 
simply, a local functional of x{t). 

We are interested in computing the probability distribution of the limiting value, = 
lim^oo x{t). In this limit, the Gaussian stochastic process, r](r), is a Gaussian variable 7700 = 77 
and the Langevin equation takes the form (to simplify notation we set Xoo = x and rj^ = if) 

V = ^ (3) 
dx 

This relation indicates that dll(x)/dx is drawn from a Gaussian distribution. We are interested 
in the distribution, p(x), of x and will try to determine it from its moments, (x l ). The partition 
function is given by the expression 

Z = / dxdr] e 2 ,5 I ^ _ _ I = / ctx — - e 2 ^ d * > (4) 

J —00 V ax j J —00 "x 

This holds only if U"(x) 7^ for all values of x. If U"(x) can vanish then the second equality does 
not hold-a well known fact of freshman calculus: we must partition the domain of integration 
into intervals, where the Jacobian is of fixed sign. The point (s) where the Jacobian vanishes are 
"excised" and the integrals computed by limiting procedure. We would like to explore in greater 
detail the symmetries of the generating function Z, written as an integral over x, when this 
limiting procedure is necessary. The reason is that this toy model will be useful as a laboratory 
for field theories. In fact it describes the boundary degree (s) of freedom of supersymmetric 
quantum mechanics in one space dimension. 
The plan of our paper is the following: 

In section [2] we show that the "classical" effective action, obtained by exponentiating the 
Jacobian, is invariant under a transformation, linear in x,whose parameter is an anticommuting 
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variable, thereby displaying worldpoint supersymmetry. Since we must introduce a regulariza- 
tion procedure, in order to integrate over the fermionic zeromodes, it is not obvious that this 
symmetry is a symmetry of the partition function, as well. We present an argument, under 
what circumstances the regularized partition function, once the cutoff is removed, will, indeed, 
respect supersymmetry. 

To gain further insight, we study the infinite set of relations between the moments, the 
zero-dimensional avatars of the Ward-Takahashi identities of field theory. 

This we do in section |3j We first establish the identities that are satisfied by the moments 
of the auxiliary field. These, however, depend quite loosely on the particular dynamical model. 
The identities between the moments of the "physical variable" are obtained by replacing the 
auxiliary field in these identities by its equation of motion. We must, therefore, show that the 
moments of the physical variable exist, for these identities to be well-defined. 

We compute these moments in two, different, ways: on the one hand we invert the "Nicolai 
map" , F = dU / cir-which can be done exactly, for the example we consider-on the other hand 
we compute the moments of the scalar variable from the moments of exp(—(dU/dx) 2 /2). 

In section H] we insert these moments in the identities satisfied by the auxiliary field and 
deduce relations for the moments of the physical variable itself. These can be evaluated either 
by using properties of special functions, or numerically. 

In section [5] we discuss the challenges that the generalization to finite dimensions poses and 
how they might be tackled, as well as the issues for describing "target space" supersymmetry 
in this formalism. 



2 Worldpoint Supersymmetry 



We would like to write 



d 2 U(x) 



dx 2 



(5) 



One way, of course, would be to write 

d 2 U{x 



dx 2 



lot 



In field theory this expression would become exp (Trlog[/"(a;)) and, since U(x) in field theory 
is a local functional of the fields, this contribution to the action would be non-local. 

On the other hand, as is well known, we may introduce it in a local way (in field theory) if 
we use two anticommuting variables, ip a , a = 1,2: 

d 2 U{x) 



dx 2 



d^d^e^^"^ 



(6) 



with 



-a/3 



-1 

1 
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The anticommuting variables introduced here are not ghosts, but as "physical" as x. There 
isn't any notion of spin so the spin-statistics theorem is vacuous. (The only way they could be 
ghosts would be if U"{x) were imaginary.) 
The effective action, therefore, becomes 

S «M) = \(^) -\^U"{x)^ (7) 

We notice a mismatch: one commuting degree of freedom versus two anticommuting degrees 
of freedom. We may restore equality by introducing another commuting variable that we will 
call, for historical reasons, F [2J: 

poo fiao 

dFe-«=i/ dFe 1 ^ 2 -^ (8) 

) J — ioo 

Here we remark that the variable F, in the field theory generalization, has an ultra-local 
propagator and that, were this propagator local, F would be a ghost. The final expression for 
the classical effective action, therefore, is 

S eS (x, ^ F) = -V + F f _ \^^, (9) 
This action is invariant under the following transformation [2j: 

5X = CaS a % 

Sip* = ( a F (10) 
6F = 

The parameter(s), ( a , a = 1,2, are anticommuting variables so we remark that this trans- 
formation mixes commuting and anticommuting variables. It is thus called a supersymmetry 
transformation. Since there is one anticommuting pair of variables, ( a , this is M = 1 supersym- 
metry. It doesn't seem to depend on the explicit form of U(x), which is the zero-dimensional 
avatar of the superpotential (this is due to the fact that we have introduced the auxiliary field, 
F, which leads to a linear realization [3j |2J). When (F) ^ 0, these relations imply that the 
fermion becomes the zero-diensional avatar of the goldstino. 

If we compute the anticommutator, generated by the transformations (flQl) . we find that 

[CQ> t)Q] = 0& {Q a , Qp} = (11) 

This relation implies that eigenstates of these operators with nonzero eigenvalue are paired. It 
is silent about the existence, or not of eigenstates of zero eigenvalue. If such latter states exist, 
then supersymmetry is realized-if not, it is broken. The resolution depends on the dynamics. 

Since it is the dynamics that will interest us, we shall work with the partition function 
and the identities of its moments and not consider the algebra itself in what follows. The 
calculations are much more direct and we need make much fewer assumptions. 
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The "equation of motion" for F, obtained by varying the effective action, is 



From eq. ([3]) we deduce that F — rj, which means that the "auxiliary variable" , F, which was 
introduced to render the supersymmetry transormations linear in x, is drawn from the same 
distribution as the noise. 

If U(x) is a quadratic function of x, then the relation between F and x is linear. This is 
the case of the free theory (provided the relation is invertible). If U(x) is not quadratic, the 
standard approach [21 El C] is to express the "physical" field, x, in terms of the "auxiliary" field, 
F through a perturbation expansion about a reference configuration. Supersymmetry implies 
identities between the correlation functions that follow from Wick's theorem applied to the 
Gaussian distribution of the auxiliary field. 

Here we would like to understand what happens to supersymmetry, in this formalism, if 
e a "U"(x) has zeromodes. 

In the next section we shall write down the stochastic identities for the case of the cubic 
superpotential, without recourse to any perturbative expansion. 

We shall use two, logically independent, methods: The first consists in solving the equation 
of motion for the auxiliary field and expressing the moments of the scalar in terms of the 
moments of the auxiliary field. 

In the second approach we shall compute the moments of the scalar from the moments of 
the classical action of the scalar. We end up with identical expressions in both cases-since we 
solved the equation of motion of the auxiliary field exactly. This is just a test case. In more 
general situations we cannot solve the equation for the auxiliary field and the second approach 
is more effective. 

We conclude with a discussion of directions for further inquiry. 

3 Worldpoint Supersymmetry beyond the classical ac- 
tion 

An example that illustrates these issues is that of the zero-dimensional model with cubic 
superpotential: 




(13) 



When we integrate out the "fermions" and the auxiliary variable, F, we find 




, 2 +\x\ 
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(In fact here there's a sign issue that doesn't have anything to do with the sign of the Jacobian, 
namely, when we integrate over two anticommuting variables, ipi and ip2 

J d^d^e^ 2 = J #!# 2 (1 + Afcfo) = -A 

This is a global sign and doesn't play any role in the calculation of correlation functions-we fix 
it once and for all.) 

We remark that U"{x) = m 2 + \x vanishes at x — — m 2 /\. If we try to compute (x p ) by 
direct sampling: 

1/2 A 

. _. f dxx p siga (m 2 + Xx) (m 2 + Xx) e~^ m x+ * x > . . 

l = _U 2 I A 

/ dx sign [m 2 + Xx) (m 2 + Xx) e 'A m X+ 2 X ) 
we see that this expression isn't well-defined, since the naive action 

S imivc = ^ (^m 2 x + ^x 2 ^j -ln|m 2 + Aa;| (15) 

is not bounded from below: In \m 2 + Xx\ — > — oo as x — > —m 2 /X. We remark, further, that the 
"classical action" , 

S = - m x H — x 

has a double well structure, with degenerate minima at x = and x = —2m 2 /X and a maximum 
at x — —m 2 /X. At this maximum the (Majorana) "mass" of the "fermion" vanishes- "chiral 
symmetry" might be realized. At the minima "chiral symmetry" is "broken" -the "fermion" is 
"massive" and its mass is equal to that of the "scalar" -supersymmetry is unbroken. If we study 
the system in perturbation theory, we can only access the vicinity of the minima and we should 
find unbroken "supersymmetry': m F = U"(x*) andm| = (U"(x*) 2 +U'(x*)U"'(x*)) = m 2 ?, since 
x* is root of U'(x*) = 0; broken "chiral symmetry" (m F ^ 0) and no hint of the instability 
at the maximum. (If we draw e~ s , we notice that it has support on the whole real axis.) 
When we try to compute the moments beyond perturbation theory, however, we will generate 
configurations around the maximum-where the Jacobian vanishes. In this case we need to 
understand, whether the theory is irretrievably sick beyond perturbation theory, or can be 
salvaged. We start to suspect that far from sickness this is a sign of health: we have tried to 
"integrate out" a "massless particle" -it shouldn't be surprising that we run into trouble, since 
we're moving along flat directions in field space. 

These are the zero dimensional avatars of the zero modes of the Dirac operator in field 
theory. So, following standard procedure, we must omit the point x = —m 2 /X, when we 
evaluate the Jacobian (i.e. when we evaluate the fermionic determinant) and "integrate over 



6 



the zeromode". Here this amounts to "excising" the point x = —m 2 /X, i.e. evaluating the 
following partition function 



Z £ = [ dx sign (m 2 + Ax) (m 2 + Ax) e h( m * x +? x3 ) + 

J-oo s v ' 

=-i 

dx sign (m 2 + Ax) (m 2 + Ax) e "^ m2x+ ^' 2 ) 

-^+e v v ' 



=+1 



and using it to compute the moments (x p ). Since we are not integrating over the entire real 
line-we have deleted the interval [— m 2 /A — e,—m 2 /\ + e]-we have broken supersymmetry 
explicitly [H However we are interested in the limit, e — > and would like to understand, 
whether, in this limit, supersymmetry is recovered, or not. 

The answer, in fact, is negative-and it isn't hard to see why. The reason is that the zero of 
the Jacobian is not, also, a zero of the action. Were this the case, then supersymmetry would 
be recovered. The proof of this statement goes as follows: 

In eq. (|T6|) the Jacobian no longer vanishes in each integral. Therefore we can perform the 
change of variables, F = m?x + (A/2)x 2 , in each and obtain the expression 



2A 



Z £ = - dFe~ F2 / 2 + dFe- F2 ' 2 = 2 dFe^ 2 ' 2 (17) 



..i^ I e^A / I e-^A 

2A ~ l ~ 2 2A ' 2 



This is the exact, regularized, partition function for the auxiliary field. 

For e finite, even if we take the limit e — > 0, with m and A fixed, the lower limit of integration 
tends to a finite value, — m 4 /(2A). This signals supersymmetry breaking, in the limit e — > 0. 
The reason is that, as e — > 0, the curve F = m 2 x + (A/2)x 2 stays fixed-and the minimum of 
the right hand side is not equal to zero. 

This, however, is easily arranged. It suffices to perform the following change of variables: 

Tfl^ A ( TTt^ \ ^ 

F = c + m 2 x+ (A/2)x 2 = c- — + - (x + —J (18) 

and fix the constant c by the condition that, in the limit e — > 0, the lower bound of the range 
of integration is at zero. In this case we immediately find that c = m 4 /(2A). For this value of 
the coefficient of the linear term, F becomes a perfect square, whose (double) zero coincides 
with the (simple) zero of the Jacobian. 

It would seem that supersymmetry has been thereby restored. Alas, this is not the case: 
having rendered the auxiliary field a perfect square means that its range is not the whole real 
axis, but, only, the non- negative part of it. Therefore the integration domain has a boundary- 
and that leads to supersymmetry breaking: 

2 f°°dFe- F2 / 2 F 
[ ' 2j co dFe- F2 / 2 * 1 } 
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We realize that the fundamental reason is that the curve, F(x), eq. (1181) . has a global maximum- 
or minimum, depending on the sign of A-that prevents us from extending the integration range 
to the full real axis. We would need it to have, at most, an inflection point. An example of this 
is provided by the quartic superpotential, 

x to? 9 A a dU n A o , . 

U(x) = — x + —x + cx <^> F = — = m x + -x 6 + c (20) 
2 4! ax 6 

If m 2 > 0, A > then the Jacobian, U"(x) = m? + (A/2)x 2 is positive definite, therefore the 
partition function becomes a Gaussian integral, centered at zero and whose integration range 
is the whole real axis. Supersymmetry is realized in the Wigner mode. 

If m 2 < 0,A > 0, then the Jacobian vanishes at two points, x = iy— 2m 2 /A. Its sign is 
positive outside of the interval (— a/— 2m 2 /A, a/— 2m 2 /A) and negative within. Therefore, the 
partition function takes the following form 



Z E = I ax I m H — x I e 2 — / ax I m H — x | e 2 + 

-oo \ 2 / ./_. /_'> m 2 V 2 

dx I m 2 H — x 2 ] e~2 J 



A /-2m 2 /AH 



- v /-2m 2 /A+ £ V ^ / ^21) 



We may perform the change of variables, F = c + m 2 x + (A/6)x 3 , in each integral. The partition 
function becomes 

•F{-y/-2m*/\-e) / .F( A /-2m 2 /A- £ ) 



oo 

oo 



F( v /-2m 2 /A+e) 



F(- v /-2m 2 /A+ £ ) (22) 



and we notice that, since F(x) is decreasing in the middle integral, F{— a/— 2m 2 /A) > F(^—2m 2 /X). 
Therefore, in the limit e — )■ 0, which is smooth, the contribution of the middle intgral cancels 
out the contribution from the overlap between the two other integrals, we are left with an 
integral over the whole real axis-and supersymmetry is restored. The linear term doesn't play 
any role at all, as far as realization or not of supersymmetry is concerned-it can only affect how 
supersymmetry is broken: For the cubic superpotential, even if its coefficient doesn't take the 
special value, m 4 / (2A), the range of integration for F still cannot cover the whole real axis and 
supersymmetry will be broken. For the quartic superpotential the linear term is completely 
invisible, as far as supersymmetry breaking is concerned. 

This analysis settles the issue of supersymmetry breaking, in principle, for these models. 
We would like, however, to understand the consequences for the moments of the scalar. 

In the remainder of this section we shall focus on the cubic superpotential and establish 
how the relations between the moments of the auxiliary field, F, even when supersymmetry 
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is broken in the way we have seen, imply, nonetheless, relations between the moments of the 
scalar. 

We start from the partition function of the auxiliary field itself, since it is now in a well- 
defined form and deduce recursion relations between its moments. These are simple enough, 
that we may solve them exactly. However they are relations between moments of the auxiliary 
field-to render them effective, we need the moments of the scalar. 

In the following subsection we establish well-defined expressions for the moments of the 
scalar variable, x, itself and use an exact inversion of the "Nicolai map", F = dll/dx, to 
express the moments of x in terms of the moments of F . The identities for the moments of the 
auxiliary field thus become identities for the moments of the scalar. However these identities 
do not carry additional information, since they simply express the Nicolai map, which is exact. 

Such information is obtained from our calculation in the next subsection 13.31 There we 
compute the moments of the scalar, without recourse to the auxiliary field at all. We obtain 
in this way expressions that involve moments with respect to an, a priori, completely different 
distribution. Nevertheless, the way we conducted the transformations implies that the moments 
thus obtained must, when substituted in the expression for the auxiliary field, lead to the 
identities that it satisfies. In our example we can see explicitly that we obtain the same expresion 
as that from inverting the Nicolai map. In more complicated cases, however, this, second, 
approach is the only one available. The simplest example is that of the quartic superpotential. 
We check our expression for the moments by substituting them in the identities satisfied by 
the auxiliary field. It is here that we really needed to compute the integrals numerically. 
We compute several identities, namely (F) = while varying the couplings and verify that 
supersymmetry is, indeed, realized, to numerical accuracy. 

The true payoff, of course, lies in the moments for the scalar, for whose distribution we do 
not have an explicit expression, but which we can, in principle, reconstruct. We can compute, in 
particular, for the cubic superpotential, the exact expression for the connected fourth moment, 
(that, in the field theory context, would control the scattering of two scalars to two scalars) 
and show that, for the cubic superpotential, it never vanishes. Therefore the distribution for 
the scalar is always non-Gaussian. 

3.1 The moments of the auxiliary field 

From the expression of the partition function for the auxiliary field, eq. (FITI) we deduce the 
following recursion relation for the moments 




) 



2 



(23) 



9 



whence we deduce that the odd moments and the even moments decouple: 



<**■■>„ = a <**-!> -b 21 

(F a ) f ,= (2(-l)(F 2 S-')> F -B 21 _ 1 (24) 

Since Z £ exists, (1) F = 1, (F) F and (F 2 ) F exist as well, all higher moments of the auxiliary 
field are uniquely specified by the first and second moments. 

These identities do not constrain the 1-point function, (F) F , whose expression is 



f 4 \ 2 



( F )f = \ . dF e~ p2 l 2 F = - " - _ F2/9 (25) 



^ 2f°° m *dFe- F2 / 

2A J C- — 



This result is interesting for several reasons: (a) A consequence of the "classical equation of 
motion", F = U'(x), for the auxiliary field, was that F = r], thus (F) = 0. This equation is, 
obviously, in contradiction with eq. (I25p . The resolution is that the symmetry is spontaneously 
broken by the interaction with the noise (parametrized by the fermions). (b) If we draw the 
distribution, 

" (F) = rfaCT (26) 

for different values of the lower limit, c — (m 4 / (2A), cf. fig. dj we notice that its maximum is at 
F* = 0, which is different from (F) F . Thus, p((F)) < p(F*). This means that, even though it 
is (F) that controls, whether supersymmetry is broken, or not, F* = is the "typical" value 
that will be drawn from p(F). The ratio, 



4^ = « p 

P(F') 



4\2 



1 



2A 



2 I 2 (°° m 4 dF e~ F2 l 2 



(27) 



We plot this ratio, as a function of the control parameter, c — (m 4 /(2A)),in fig. [2j The result 
is that, for c - m 4 /(2A) w -0.3, the ratio falls to ~ 0.75. For |c - (m 4 /(2A))| > 1, this ratio 
is, practically, equal to one. But, for these values of the control parameter, p(F) is, also, very 
small, so large samples are required. 

Therefore, if c = m 4 /(2A), in which case = 0, we would be tempted to conclude that 
supersymmetry is realized, even though it is spontaneously broken, if we couldn't distinguish the 
typical from the average value. On the other hand, this might, also, imply that supersymmetry 
breaking effects are much harder to detect than expected, due to this background. This issue 
certainly deserves further study in more realistic cases. In disordered systems this was noted 
in ref. [H]; in ref. [15] it was shown that it is, indeed, possible to distinguish the typical from 
the average value of certain correlators in disordered spin chains. If c ^ m 4 /(2A), then E>k ^ 
and the identities satsfied by the moments of the auxiliary field acquire "anomalous" terms. 
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Figure 1: Some examples for the distribution of the auxiliary field, p{F), vs. F > (c — 
(m 4 /(2A)), for different values of c - (m 4 /(2A)). (a) c - m 4 /(2A) = -10-here F* = « (F); 
(b) c - (m 4 /(2A)) = -1-here F* = ^ (F), F* is the typical value, (c) c = m 4 /(2A)-here, 
too, F* = ^ (F). (d) c — (m 4 /(2A)) = 1-here F* = doesn't lie in the sampling domain, 
F* ^ (F) ^ 0. 




Figure 2: p((F))/p(F*) as a function of c - (m 4 /(2A)). 
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Nevertheless, from the relation between the auxiliary field and the scalar, eqs. flH]) lead to 
constraints on the moments of the scalar. 

We shall present two approaches for computing its moments, (x p ) , none of which entails an 
expansion in any coupling constant, but is exact, (to machine precision). The first simply relies 
on the fact that we can solve quadratic algebraic equations explicitly; the second dispenses with 
this requirement, but relies on the ability of evaluating integrals numerically. 



3.2 Inverting the Nicolai map 

The relations obtained above are quite model independent-the existence of the boundary limit 
of integration, B, doesn't depend on a specific model, only on a class as a whole-to give them 
content we must replace the auxiliary field by the derivative of the superpotential under study 
and define how we compute the moments of the scalar. 

In this subsection we shall assume that we can solve the equation of motion of the auxiliary 
field, 

p= dU_ 

dx 

exactly. This is what is meant by "inverting the Nicolai map" . We shall use this knowledge to 
write the regularized moments, (x p ) £ of the scalar in a form that is particularly suited towards 
numerical calculation and shall try to check, whether these moments, when inserted in the 
expression ((dll/dx) q ) £ , in the limit e — > 0, lead to identities consistent with Wick's theorem. 

Since we are particularly interested in the role of the linear term of the superpotential, we 
shall include it, with an arbitrary coefficient. 

The expression for the regularized moments of the scalar variable, x reads as follows: 

- r^~ £ dx (m 2 + \x)xPe-^ c+m2x+ ^ x2 ¥ + r m 2^ dx (m 2 + Xx)x p e-^ c+m2x+ ^ x2 T 

I n \ — OO J — — \-£ 

(x p ) £ = - 



jZJ~ £ dx (m 2 + Xx)e-^ c+m2x+ ^ x2 T + JZ* dx (m 2 + Ax)e~K c + m2 *+^ 2 ) 



We perform the change of variables, 



(2f 



F = c + m 2 x + —x 2 dF = dx(m 2 + Ax) 



whereupon the moments are computed as 



P° m 4 2A dFe~ F2 ' 2 [x-(F¥ + x+(F) 

( x p\ = Jc ~^ + — . - ( 2 Q) 

{ h 2JZ^ +s ^dFe- F2 / 2 [ ] 



where 



x ± (F) 




m 4 \ 



F - {C ~2X) 



(30) 
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We end up with the following compact expression for the regularized moments, (x p ) £ : 



2\p p / \ f°° 4 2, dFe~ F2/2 (f - c+ ~ '' N V '' " 

(*">« = ( -4 ) 2 £ ' p ^ Jc "* +¥ [m { 



' fc=0,fceven N ' "2A"" 1-- 2~~ v^-U 

?)Vr£(> )(§)>), 

where (F k ^ are computed using the distribution ff26]) . 

From these expressions we can, in principle, recover the density, p(x) and the action, S(x). 

It must be stressed that (x l \ are the exact moments. They are the coefficients in the Taylor 
series expansion of the Fourier transform, p(k), of the density, 

P -S{x) 

p{x) 



J — oo 



around k = 0: 

p(k) = I dxe lkx p(x) <^> i 



iA- >■ ; ^ / dl p( k ) 



dxx 1 p{x) = (x') 



r\h l 

. i k=0 (32) 

1=0 ' J ^ 00 



While the (exact) moments in eq. ( 13TT) are certainly well defined quantities, whether they do, 
indeed, define a density, p(x) and, thus, an "action", S(x) = — lnp(x), for all values of the 
"coupling", 2\/\m 2 \ 2 , is a question that deserves further study, that will be reported elsewhere. 

In these calculations we relied on the fact that we could solve F = dU/dx exactly. There 
are, however, situations when we cannot invert the Nicolai map as efficiently as done in this 
case. An example would be that of a quartic superpotential or many superfields. Then we 
would have to solve a cubic equation, or a system of non-linear equations-and this cannot be 
done as easily as for one quadratic equation. So we need to develop methods for computing 
the moments of the scalar that do not rely on the Nicolai map. 

To test the approach we shall apply it to the case at hand, where a direct comparison is 
possible and to the quartic superpotential. 



3.3 A new measure for the scalar 's moments 

Our calculations were based on inverting the "Nicolai map" , u = dU / cfe-however this inversion 
didn't involve any approximations at all, since we were "lucky" that we had to solve a quadratic 
equation, whose roots are known analytically. Had we chosen not to invert it-or could not invert 
it exactly-we would have to compute the moments, (x p ) and use them to compute the moments, 
(F q ), of the auxiliary field. Of course we could have computed the moments in a perturbative 
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expansion [21 [7], but this would weaken considerably the conclusions we could reach. With 
numerical methods available we can explore what happens beyond perturbation theory. 

It is instructive to compare the two approaches in this case. 

The idea is to compute the moments (x p ) in terms of the moments of 



P(y) = exp 
Our starting point is eq. (|28|) . We set 

and write the numerator as follows: 



A , 

C "2A + 2^ 



-S(y 2 ) 



(33) 



y = x + 



T 



A 



dy e 



S(y 2 ) 



ni 



2\P 



y y 



A 



+ 



dye 



S(y 2 ) 



m 



2\P 



y y 



A 



(34) 



The denominator is obtained setting p = in this expression. Once more, the e — > limit is 
smooth, so the moments (x p ) may be computed by 



lim (x p ) £ = lim 



JZ dy e- s ^y (y-^) P + r dy e~ s ^y (y - ^ 



£^0 

S~dye< c -^)\ 



(35) 



J™ dye 



The numerator can be simplified: 

,2\P 



-y 



m 

T 



y 



2 - 



m 



2 \ P 



A 



m 

T 



2\P P 



E 

k=0 



E 

fc=0 



V 
2k 



A 



m 

T 

2k 



2 \ P~ k 



[(-y) k + y k ] 



(36) 



y 



2k 



This leads to the change of variables 2ydy = d(y 2 ) in the numerator. After some further 
straightforward algebra, we end up with the expression 



2\P 



fc=0 



p 

2k 



2X\ k f~dye 
Jo d V ( 



(37) 



which, we perceive, after the fact, is identical to eq. (13T|) ! We didn't need to solve the equation 
F = dU/dx—we only used the knowledge of the root, x = — m 2 /A, of d 2 U/dx 2 = m 2 + Ax = 0, 
around which we needed to impose the excision. 
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For the quartic superpotential, with m 2 < 0, A > 0, the procedure is the following: We 
write U"(x) = m 2 + (A/2)x 2 = (A/2) (a; — xi)(x — x 2 ), with x\ < x 2 . If we define the "reduced 
partition function", ( a> b, by the expression 



c a , 6 = r<ke-5(^ s 

J a 



(38) 



the moments, (x p ), are given by the expression 

(^-oo.oc 



(x*> 



IU") 

C-oo,xi-e (U x l -oo,xi-e ~~ Cxx+e,x 2 -e (U xP ) Xl +s,x 2 -e Cx 2 +e,oo (U xP ) X2 +e,oo 



lim {x^i — lim 

£-5>0 e->-0 C-oo,xi-e \U") _ OO Xl _ £ — C, Xl +e,x 2 -e (U") Xl+£ ^ X2 - e + Cr 2 +£,oo (^")x2+£,oo 

(39) 

where (C) a6 denotes average with respect to Co,6- Of course here this is just a convenient 
shorthand-it is more efficient to use Simpson's rule than Monte Carlo (we've used both). But 
for the more general cases, these expressions generalize better. They are "templates": for 
concrete applications it is useful to try to simplify them, in order to generate efficient code, as 
was, already, done, for the cubic superpotential. 

In the following section we shall study the identities that these moments satisfy. 



4 Moments and Identities 

Let us start with the cubic superpotential. 

The expressions for the moments display, first of all, scaling: 

<MxN^) (4o) 

The prefactor, (m 2 /\) p indicates the "canonical dimension" of the p— th moment-which co- 
incides with its "scaling dimension". The "scaling function", / p (-), is a function only of the 
combination, A/m 4 . The combination m 2 /A simply sets the scale, it is the combination A/m 4 
that plays the role of the coupling constant. In other words, the combinations that are relevant 
for describing the moments are not m 2 and A, but m 2 /A and A/m 4 . 

The moments are perfectly regular functions of the cutoff, e, and have a smooth limit, as 
it is removed, e — > 0. The coefficient, c, of the linear term of the superpotential simply affects 
the value of integration limit-even if this can be pushed to infinity in the denominator (thanks 
to the factor of 2), it survives in the numerator. However, if 
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then the identities satisfied by the auxiliary field, F imply that supersymmetry is spontaneously 
broken: (F) F 7^ but all other identities that express the Gaussian distribution for F are 
exactly satisfied. 

If c 7^ m 4 / (2A), then the fermion is still the goldstino, but the WT identities have anomalies. 
If we compute (F), in the limit e — » and for c = m 4 / (2A), we find that 



,„ m 4 .,. A,,, m' m' m' J^dFe- F '/ 2 F 
< F > = 2A + "' <*> + 2 <^> = . 2A--+2A + 2 FdFr"» 



(42) 



=0 

7^0 



identical with the expression (1251) . when c = m 4 /(2A). 

Similarly, if we compute (F 2 ),when c = m 4 /(2A), in the limit e — >■ 0, we find 



A 2 / m 2X 



^((.V4(,^ + 6 (r>(^) 2 + 4(.)(^) 3 + (^) 4 



(43) 



Substituting the expressions for the moments ( 13T|) . we find, indeed, that (F 2 ) = (F 2 ), i.e. that 
the right hand side of eq. (1151) is equal to 1. 

These are but two examples of eq. (123|) . which implies that the moments, (x p ), of the scalar, 
given by eq. (15T1) or (155]) . satisfy the following identities 

where B k is defined by eq. (j2"5|) . These identities would be quite difficult to guess, if the only 
inputs were the moemnts, (x p ), given by eqs. ( 15T|) or (137|) . 

The stochastic approach, therefore, helps us (a) to obtain concrete expressions for the 
moments, (x p ), of the "physical" variable, x, without prior knowledge of the probability density, 
p(x), and (b) to deduce identities between these moments, that would be very difficult to deduce 
from scratch. 

An interesting question is, whether the moments are compatible with a Gaussian distribution 
for the scalar. One way to address it is to compute the connected four-point function, 

<x 4 ) c = (x 4 ) - 3 (x 2 ) 2 + 3 (x 2 ) {x} 2 - (x 3 ) (x) (45) 

which vanishes for a Gaussian distribution as a consequence of Wick's theorem. Using eqs. (l3T]) .or 
( )35|) . we find, indeed, that it doesn't vanish. By dimensional analysis and scaling it is of the 
form 
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and we expect that f(g) — > as g — >■ 0, since, for A/m 4 — )• the superpotential is quadratic. 
Using either of the expressions (I3ll . (l35l) we find that 

for the case when c = m 4 /(2A). We can, of course, compute it for any value of the coefficient of 
the linear term. We conclude that, for any, finite, value of A/m 4 , the distribution of the scalar 
differs markedly from that of a Gaussian, and the "non-Gaussian" effects become more marked 
with increasing A/m 4 at fixed m 2 /A. 

Now let us turn to the quartic superpotential, when m 2 < 0, A > 0. As we saw, despite the 
presence of fermionic zeromodes, supersymmetry is realized. However since we can't invert, 
efficiently, the Nicolai map, but compute the moments of the scalar by a way that doesn't, 
immediately, seem related to the auxiliary field, to show that the moments, computed this way, 
do satisfy the stochastic identities is not completely obvious. 

The master identity, that expresses the fact that the auxiliary field, F = m?x + (A/6)x 3 , in 
the limit when the cutoff is removed, is drawn from a Gaussian distribution, reads 



k{F k ~ l ) F = {F k+1 ) 



F 



k I — |m \x H — x ) = ( ( — |m \x H — x 

6 ; / \v 1 1 e 

along with the "initial condition", (F) = 0. The master identity, for k — 1, implies the first, 
non-trivial, relation 



<**)„= <!>„ = !=►/( 




(49) 



We shall start with the 1-point function for the auxiliary field: Supersymmetry implies that 
(F) F = 0, therefore, 



A 

\m , 

6 



2| (x) + -(x 3 ) = (50) 



Once more it's useful to work with rescaled variables. Let us set 

01 21 \ 1/2 

2 77T 1 X 



xsX { -ir ) (51) 

The moments then are calculated through 

•2|m 2 |V /2 l Jvf i , v2 ^ vv -^(-x + 4) 2 _ /2|m 2 |V /2 ((-1+X 2 )X^) 

((-1+^ 2 )) 

(52) 



{XP) = (^f) \ j°°JX{-l + X 2 )XVe 
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The scaling relation for the moments of the quartic superpotential becomes 

and we can write the master identity (|48l) in terms of the rescaled variable X, in the following 

way 

X 3 \ k+1 \ ( \ \ 2 If X 3 \ kl \ I ( X 3 \ 2 \ / / X 3 ^ k 



(54) 



with 



3 / / ™ |m2|3 ( v , X 3y 

° 7 ' /_°^rfX(-l + X 2 )e-^V- x+ -J 



The relation, = becomes 



3 

If we write it out: 



-(X) + ^p- = (-X + ^-\ = Q (■-)()) 



dX -X + — (-1 + X 2 ) e ~ J = (57) 



it becomes obvious, since, in this guise, the integral exists, thanks to the exponential factor 
and the integrand is the product of an even function and an odd function, integrated along an 
interval, symmetric about the origin, therefore it vanishes. Indeed, any odd moment vanishes 
as a result of the symmetry, X -h- —X. 

Let us look at a, possibly, less trivial relation, namely {F 2 ) F = 1. In terms of the (X p ), it 
reads 

/ ,\2 |m 2 | 3 I Vl x 3 \ 2 



3 / / ?\m 2 \ 3 ~, I™ 2 ! 3 ( v^x3\ 2 2lw 2 l 3 

(58) 

If we could perform the change of variables, F = \J\m 2 \ 3 /\ (— X + (X 3 /3)), it would be obvious. 
The analysis in the previous section confirms this expectation. Let us show that it is possible 
to recover this result by computing the moments according to eq. f|52|) . 

If we draw exp ^— (|m 2 | 3 /A) (—X + (X 3 /3)) 2 j for different values of the coupling constant, 

A/|m 2 | 3 , cf. fig. |3] that cover the interval from "weak" (A/|m 2 | 3 < 1) to "strong" (A/|m 2 | 3 > 1) 
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Figure 3: exp ( — (|m 2 | 3 /A) (—X + (X 3 /3)) 2 J as a function of X for different values of |m 2 | 3 /A. 
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coupling, we note that, for weak coupling, the density is concentrated around the "classical" 
minima, while, as the coupling increases, it becomes easier to explore other values, in particular 
around the local maxima of the "classical action" -however we do note that, for \X\ > 5, the 
density becomes negligible: for |X| » 5 the density is smaller than ~ exp(— 37 2 |m 2 | 3 /A). 
Therefore we may cutoff the integrals at |X| = 5. 

In fig. H] we display, for different values of A/|m 2 | 3 , the absolute value of the difference 
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Figure 4: Absolute value of the difference of the two sides of relation fl58l) as a function of 
/A. The vertical scale is logarithmic. The different curves correspond to different values 
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of points for Simpson's rule, used for evaluating the integrals, in the interval —5 < X < 5. 



between the left hand side and the right hand side of eq. fl58l) . We find, once more, a result 
that is consistent with zero, to machine precision, when numerical integration errors become 
negligible. Our code passes the test. 
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5 Conclusions and Perspectives 



We have studied the partition function of probabilitiy distributions of one variable, that are 
generated by solutions of Langevin's equation. This defines a zero-dimensional field theory. 
We have shown that the induced probability distribution enjoys worldpoint supersymmetry, 
that can be broken by surface terms, induced by fermionic zeromodes. The "fermions" enter 
through the noise term. 

While a linear term in the superpotential, also, contributes surface terms, in fact it cannot 
cancel their contribution. In this way we provide concrete examples for the boundary terms 
studied by Parisi and Sourlas [3J. 

Nevertheless, this breaking doesn't imply that supersymmetry isn't useful: as we have 
shown, even broken, it, still imposes relations on the moments and all higher moments of the 
auxiliary field are determined by its first and second moments. So the auxiliary field is still 
drawn from a Gaussian distribution, "deformed" in an interesting way, due to the interactions. 

For practical purposes, we can find examples, where supersymmetry is spontaneously bro- 
ken, but the goldstino is "hard to find" , since the average value of the auxiliary field is not the 
"typical" value, that is drawn from its distribution. 

The scalar field, on the other hand, isn't drawn from a Gaussian distribution-but we can 
compute its moments, at any value of the coupling, quite explicitly and they determine its 
distribution (its effective action, in the field theory context), in principle. We have determined 
exact expressions for the moments of the cubic superpotential by two independent methods: 
inverting the Nicolai map and by establishing the relation to the moments of the scalar part 
of the classical action. It is this second approach that most naturally generalizes to higher 
dimensions, where numerical simulations are the natural tool, which we tested on the quartic 
superpotential, where supersymmetry is realized, whatever the sign of the quadratic term of 
the superpotential. 

So the situation is by no means "all or nothing": supersymmetry either is realized and 
constrains everything or is broken and is not useful, but is much more interesting. The in- 
formation that has not been appreciated lies in the identities between moments or correlation 
functions. While in the "traditional" approach to supersymmetric field theory, one starts from 
the "physical" (albeit classical) fields and the auxiliary fields are, indeed, quite hard to find [7], 
in the stochastic approach the auxiliary fields are the starting point and the "physical" fields 
are "emergent" quantities. In this approach supersymmetry is a natural symmetry. So to un- 
derstand it better, it might seem that the stochastic approach, through the access it provides 
to the moments/ correlation functions, might prove much more useful. 

In fact a natural direction of further study is the following: The reason why the auxiliary 
field is drawn from a Gaussian distribution is that it was, initially, at least, identified with the 
noise in the Langevin equation. So we might consider other stochastic equations, where the 
noise is drawn from other distributions, e.g. "colored" noise [12] and, of course we can imagine 
further generalizations. These might lead to more general supersymmetry-breaking patterns, 
that definitely deserve to be studied in higher dimensions. 
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The most straightforward generalization is to "worldline supersymmetry" , i.e. supersym- 
metric quantum mechanics [9j [10J ITT] . This has been, mainly, studied in a way that stresses 
its similarity to (non-relativistic) quantum mechanics rather than quantum field theory (an 
exception is ref . [10] , where, however, the focus was on topology and not on the stochastic iden- 
tities themselves, that were studied in ref. [7J in an approach that combined canonical and path 
integral approaches ). Lattice studies [J7] use the Nicolai map but, apparently, do not consider 
the stochastic identities themselves. These are "broken" by the lattice since the auxiliary field 
has a local propagator, that becomes ultra-local only in the continuum limit (at least for the 
free theory). 

In this paper we have tried, on the contrary, to set up the formalism that most readily 
generalizes to quantum field theory, namely the (zero-dimensional analog of the) path integral 
formalism and the WT identities between moments. We have thereby elucidated the role of 
the distribution of the auxiliary field in characterizing supersymmetry breaking and developed 
the tools that will be useful for the higher dimensional cases. 

Our results can be considered as describing the degrees of freedom that live on the boundary 
of Euclidean time, with supersymmetric quantum mechanics describing the "bulk" theory. For 
the cubic superpotential, supersymmetry is therefore broken both in the bulk 0, El El El II] 
and on the boundaries. For the quartic superpotential it is broken only in the bulk. These 
statements, however, do not take into account the coupling between bulk and boundaries, which 
is the issue that remains to be studied (work in this direction, in the canonical formalism, was 
done in ref. [TT]). 

We have studied the case of one chiral multiplet and a natural generalization is to Nf > 1 
chiral multiplets-which also lead naturally to theories with extended supersymmetry. In that 
case the technical difficulty stems from the fact that we must "excise" not just a point but 
a vector space, generated by the zeromode(s) of the Jacobian, Nf x Nf, matrix. The use of 
collective coordinates [TB] is the natural tool, already at this level. 

Finally, another topic concerns supersymmetric gauge theories. Here the problem is how 
to handle gauge invarince. In stochastic quantization [18] gauge fixing is not mandatory, but 
only gauge invariant quantities are well-defined when the equilibrium limit is taken. If we 
work on the lattice, with gauge links that take values in a compact group, all observables are 
gauge invariant, but the gauge links are constrained: e.g. they are unitary matrices, for unitary 
groups. We may solve these constrants, however [TjJ], so it will be interesting to study the 
corresponding Langevin equation and whether we may obtain supersymmetric gauge theories 
this way. 

For, indeed, the question remains, how "target space" supersymmetry can be obtained 
by the stochastic approach. Parisi and Sourlas did succeed in obtaining the Af = 2 Wess- 
Zumino model in two spacetime dimensions, but found an obstruction in four dimensions, de 
Alfaro, Fubini , Furlan and Veneziano [7J obtained some, partial, results, in particular for gauge 
theories, but the issue remains open, how to write a lattice theory in this formulation. This is 
one reason why the "orbifold" approach [20] has since been used. It is, therefore, of interest to 
test its assumptions by another formulation. 
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